Critical behavior and scaling in trapped systems 
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We study the scaling properties of critical particle systems confined by a potential. Using 
renormalization-group arguments, we show that their critical behavior can be cast in the form 
of a trap-size scaling, resembling finite-size scaling theory, with a nontrivial trap critical exponent 0, 
which describes how the correlation length £ scales with the trap size I, i.e., £ ~ 1° at T c . 8 depends 
on the universality class of the transition, the power law of the confining potential, and on the way 
it is coupled to the critical modes. We present numerical results for two-dimensional lattice gas 
(Ising) models with various types of harmonic traps, which support the trap-size scaling scenario. 

PACS numbers: 05.70.Jk, 67.85.-d, 64.60.F-. 



Recent experimental developments, namely the 
achievement of Bose-Einstein condensation in dilute 
atomic vapors of 87 Rb and 23 Na [H, 0], have attracted 
great interest in the study of atomic systems in a trap- 
ping potential. In the Bose-Einstein condensation sce- 
nario, a macroscopic number of trapped bosonic atoms 
accumulate in a single quantum state and can be de- 
scribed by a condensate wave function, which naturally 
provides the order parameter of the phase transition. The 
critical behavior of a trapped Bose gas has been recently 
investigated experimentally [3j, observing an increasing 
correlation length compatible with the behavior expected 
at a continuous transition in the three-dimensional XY 
universality class, see, e.g., Ref. [H and references therein, 
which is also the universality class of the superfluid tran- 
sition in 4 He, see, e.g., Ref. However, the inhomo- 
geneity due to the trapping potential strongly affects the 
phenomenology of phase transitions observed in the ab- 
sence of a trap. For example, correlation functions of the 
critical modes are not expected to develop a diverging 
length scale in a trap. Therefore, a theoretical descrip- 
tion of the critical correlations in systems subject to con- 
fining potentials, and of how they unfold approaching 
the transition point, is of great importance for experi- 
mental investigations of the critical behavior of systems 
of trapped interacting particles. 

We consider a trapping potential 



U(r) 



(1) 



where v and p are positive constants and I — w _1 is the 
trap size, coupled to the particle number. Harmonic po- 
tentials, i.e., p = 2, are usually realized in experiments. 
The effect of the trapping potential is to effectively vary 
the local value of the chemical potential, so that the par- 
ticles cannot run away. Let us consider the case in which 
the system parameters, such as temperature, pressure 
and chemical potential, are tuned to values corresponding 
to the critical regime of the unconfined system, where the 
correlation length diverges as £ ~ t~ v and the correlation 
function behaves as G(r) ~ \/ r d-2+v at Tq (p = T/T c -1 
is the reduced temperature and v and n are the critical 



exponents of the universality class of the transition) . In 
the presence of a confining potential, the critical behavior 
of the unconfined homogeneous system could be observed 
around the middle of the trap only in a window where £ 
is much smaller than the trap size but sufficiently large 
to show the universal scaling behavior. If £ is large but 
not much smaller than the trap size, the critical behav- 
ior gets somehow distorted by the trap, although it may 
give rise to universal effects controlled by the universal- 
ity class of the phase transition of the unconfined system, 
similarly to finite-size scaling in homogeneous systems of 
finite size [H, @] . In the present paper we investigate the 
critical behavior of trapped systems, putting on a quanti- 
tative ground the above qualitative scenario, and in par- 
ticular the relevant effects of the confining potential. Us- 
ing renormalization-group (RG) arguments, we show that 
the critical behavior of the trapped system can be cast in 
the form of a critical trap-size scaling, resembling stan- 
dard finite-size scaling theory for homogeneous systems 
at a continuous transition [y, |2j , but characterized by a 
further nontrivial trap critical exponent. 

For the sake of demonstration, as a simple model of 
trapped particles, one may consider the d-dimensional 
lattice gas model defined by the Hamiltonian 



W L gas = X! Pip i ~ p X! Pi + 2U ( r i)P* 

(ij) i 1 



(2) 



where the first sum runs over the nearest-neighbor sites 
of the lattice, pi = 0, 1 depending if the site is empty 
or occupied by the particles, p, is the chemical poten- 
tial, and U(r) is the potential ([I]). Far from the origin 
the potential diverges, thus the expectation value of the 
particle number tend to vanish, and therefore the parti- 
cles are trapped. The lattice gas model ^ can be ex- 
actly mapped to a standard Ising model, by replacing 
Si = 1 — 2pi, thus Si = ±1, obtaining 



H = - J SjSj - h Sj - }^ U(n)sj 
(ij) % % 



(3) 



where h — 2qJ + p/2 (q is the lattice coordination num- 
ber). In the absence of the trap, this model shows a criti- 



2 



cal behavior characterized by a diverging length scale, at 
the critical point T = T c and h = h c = 0. Critical cor- 
relations do not develop a diverging length scale in the 
presence of the confining potential, i.e., at fixed v > 0. 
We want to study how the critical behavior is distorted 
by the trap, and how it is recovered in the limit v —> 0. 

Our starting point is a scaling Ansatz which extends 
the scaling law of the RG theory of critical phenomena Q 
(see also Refs. 0, HI), to allow for the confining potential 
(fT|). We consider a standard general scenario in which 
the transition of the unconfined d-dimensional system is 
characterized by two relevant parameters, t = T/T c — 1 
and h, as in models (J2E]). We write the scaling law of 
the singular part of the free energy density as 

F(u u u h ,u v ,x) =b- d F(u t b y \ u h b yh ,u v b y »,x/b) (4) 

where b is any positive number, x is the distance from 
the middle of the trap, Utj t . v are scaling fields associated 
with t, h, and v respectively. They are analytic functions 
of the system parameters, behaving as u t ~ t, Uh ~ h, 
and u v ~ v when t, h,« — ► 0. yt.h,v are the corresponding 
RG dimensions: yt = 1/v and ijh = (d + 2 — rj)/2, while 
y v must be determined (see below). We are neglecting 
irrelevant scaling fields, because they do not affect the 
asymptotic behaviors. Then, fixing u v b Vv — 1 and intro- 
ducing the trap size I = v , we obtain 



F = r 0d F(utl eyt ,u h l 0Vh ,xr e ) 



(5) 



where 9 = l/y v is the trap exponent. From Eq. ([5|) one 
can derive the trap-size scaling of other observables. A 
generic quantity S is expected to behave asymptotically, 
when \t\ — > and in the large-Z limit, as 



s = r 6y °f s {ti e ' v ) = r ey °J s {£,i- e ), 



(6) 



where y s is its RG dimension, f s and f s are universal 
functions (apart from normalizations). £ is the correla- 
tion length, which behaves as £ ~ I at T c , and as £ ~ t~ v 
when I — > oo. Note that f s (x) ~ x~ y " for x — > 0, so that 
5 ^ when Z — > oo, which is the scaling behavior in 
the absence of the trap. These results are quite general, 
they apply to the lattice gas model ([2]), as well as to more 
general interacting gas systems, fluids, etc... 

The trap exponent 9 can be computed by analyzing 
the RG properties of the corresponding perturbation at 
the critical point. In the case of the model (jl]), it can be 
represented by Py = J d d xU(x)4>(x), where 4>(x) is the 
order-parameter field of the (/> 4 theory which describes the 
behavior of the critical modes, see, e.g., Ref. Since 
the RG dimensions of the potential U (x) and the field <fi, 
respectively yxj = py v - p and y (j> =(3/v = {d-2 + rj)/2, 
are related by yu + y<p = d, we obtain 



1 



2p 



Notice that 9 — > 1 when p — > oo, which is consistent with 
the fact that when p — » oo the effect of the trapping po- 
tential is equivalent to confining a homogeneous system 
in a box of size L = I with fixed boundary conditions, 
thus leading to a standard finite size scaling where the 
RG dimension of the size L is formally minus one [H, 0] • 
As examples of observables in model ([3]), we consider 
the local magnetization and the energy density in the 
middle of the trap, i.e., 



M = (s ), 



E = (s Sl 



(8) 



(where 1 is a nearest neighbor of the center site 0), and 
the correlation function 



Go(r) 



{s s r ) 



(So)(Sr) 



(9) 



d + 2 - r] + 2p 



(7) 



In the lattice gas model M and Go (r) are related to the 
particle density and its correlations. Their asymptotic 
trap-size scaling behaviors are given by 

Mo = i- ep/v f m {ti 9/lJ ) = r e P /u fm{y- e ), (io) 
G (r) = i- e y g {ti e '\ri- 9 ) = r e "/ ff (£rV/£), (ii) 

E = E ns (t) + l-^ d - 1 /^f e (tl e ^), (12) 

where 0/u = (d — 2 + rj)/2 and E ns (t) is a nonsingular 
function. The scaling behavior at T c can be obtained 
from the above equations by setting t = 0, 

M ~ r fl ^, Eo-Ecrvl-W- 1 /"), (13) 
G (r) ~ r%(rZ~ e ). (14) 

In order to check the trap-size scaling scenario, we con- 
sider the square-lattice Ising model J3|) with J = 1 (for 
which we know the critical temperature and energy, i.e., 
T c = 2/ln(y/2 + 1) and E c = V2, and the critical expo- 
nents, v — 1, r\ = 1/4) in a harmonic trap, i.e., the poten- 
tial (fT]) withp = 2. According to Eq. 0, the correspond- 
ing trap exponent is 9 = 16/31. Beside two-dimensional 
(2D) traps, we also consider one-dimensional (ID) traps, 
where the confining potential acts only along one direc- 
tion, i.e., U(r) — v p \x\ p depends only on x, while there is 
translation invariance along the other y direction. Note 
that the trap-size scaling formulae, including the value of 
the trap exponent 9, apply to both 2D and ID traps. 

We performed several Monte Carlo (MC) simulations 
around the critical point and for various values of the trap 
size I. The lattice size L and harmonic potential of the 
simulated systems were chosen to have the spin variables 
effectively frozen at the boundary, making unnecessary 
the use of larger lattices [l(| • The MC results for Mq and 
Eq at T c provide accurate checks of the trap-size scaling 
predicted by Eqs. (US]). In Figs. [T] we plot M = l 2/31 M 
and AE = I (Eq — E c ) versus l^ 1 , which is the order of 
the expected leading corrections to the asymptotic scal- 
ing ll|. The lines represent fits of all available data, i.e., 
for L — 2™ with n — 3, ...,8, to Mq = a + b/l (in all cases 
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FIG. 1: Trap-size scaling of Mo and Eo at T c : we plot Mo = 
l 9r > /2 M (below) and ~AE = l e (E -E c ) (above) vs r\ for the 
model ^ with ID and 2D harmonic traps (denoted by "a") 
for which 6 = 16/31, and the model (|15[) with a 2D harmonic 
trap (by "6") for which 6 = 2/3. The dotted lines show linear 
fits. 
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FIG. 2: Trap-size scaling of Mo for the model (J2]| in a 2D 
trap. Mo = l 2/31 M and L = ^2L 



with x 2 /dof < 1). [12j In Fig. [2] we check the scaling of 
Mq around T c in the case of the 2D trap, by plotting 
M = l 2/31 M versus tl 16 / 31 . As predicted by Eq. flQ|). 
the data approach a scaling function / rn (^ 16 ^ 31 ) in the 
large-^ limit, where scaling corrections are suppressed. 
Fig. [3] shows results for the correlation function ([9]) in the 
case of the 2D trap, at T c . They clearly support Eq. (fl"4|) . 
indeed the data of G (r) = l 4 / 31 G (r) follow the same 
scaling function g(rl~ 16 / 31 ) for all trap sizes (scaling cor- 
rections are very small and not visible in Fig. [3]). Go(r) 
appears to decay rapidly at large distance, with a length 
scale behaving as £ ~ I 6 . 

To further check the trap-size scaling scenario, we also 
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FIG. 3: Scaling at T c of Go(r) for the model © in a 2D trap. 
G (x) = l 4/31 G (x) and L = V21. 



consider another Ising-like model with the confining po- 
tential (Q} coupled to the energy density, i.e., 



(ij) 



U(rij)]i 



(15) 



where the first sum runs over nearest-neighbor sites, 
and Tij = (r, + rj)/2. Also in this case the spin 
variables get frozen at large distances. More precisely, 
lim^_,.o+li m |r|-too( s r) = 1 at any T, corresponding to 
vanishing particle density in the lattice gas model. The 
effects of the confining potential at the transition can be 
again described by a trap-size scaling, cf. Eqs. (|10M14|) . 
but with a different trap exponent 9. Indeed, a RG anal- 
ysis of the perturbation arising from the potential U in 
model (p]l. On = J d d x U{x)(f> 2 (x), gives [3 



pv 



l+pv 



(16) 



We again consider a square-lattice model with a harmonic 
2D trap. Since v = 1, Eq. |QU) gives 9 = 2/3, which 
differs from the value 16/31 for the model ([3]). We per- 
formed MC simulations similarly to the previous case. 
Results for the trap-size scaling of the magnetization Mq 
and the energy Eq at the origin, and the correlation func- 
tion ([5]) are shown in Figs. [TJ 2] and [5] They again fully 
support the predicted trap-size scaling behaviors. 

In conclusion, we have shown that the critical behavior 
of systems in a confining potential can be described by a 
universal trap-size scaling (expected to be largely inde- 
pendent of the microscopic details of the model), char- 
acterized by a trap exponent 9 which describes how the 
correlation length scales with the trap size, £ ~ l e at T c . 
The exponent 9 essentially depends on the universality 
class of the transition, the power law of the confining po- 
tential, and on the way it is coupled to the critical modes. 
These results are very general, and by no means limited 
to two dimensions. 

We finally discuss the critical behavior of a 3<i in- 
teracting Bose gas trapped by a harmonic potential, 
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FIG. 4: Trap-size scaling of Mo for the model ([To)) in a 2D 
trap. M = l 1/12 M and L = ^2/31. 
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FIG. 5: Scaling at T c of G (r) for the model (JT5J) . G (a;) 
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which has been recently investigated experimentally 
This system is expected to undergo a continuous tran- 
sition in the 3c? XY universality class, characterized by 
a complex order parameter <p(x) with U(l) symmetry. 
The confining potential U{x) is coupled to the particle 
density. The corresponding perturbation is [HI Qu — 
J d 3 x (\x\/l) 2 \(j)(x)\ 2 , where 0(x) is the order-parameter 
field. The same RG arguments which led us to Eq. (fl"6"l) 
give 9 = 2u/(l + 2v), thus 9 = 0.57327(4) using |l| 
v = 0.6717(1). This implies that one can neglect the 
trap effects only when the correlation length satisfies 
£ < I 9 = I - 573 . The experimental results of Ref. 
were likely taken in this region, and led to the estimate 
v = 0.67(13), by fitting the data to the standard behav- 
ior £ ~ t~ u . For larger values of £, when £ and l e become 
comparable, the trap-size scaling discussed in this paper 
is expected to provide the correct critical behavior. We 
believe that experiments probing the trap-scaling regime 
can be very interesting, analogously to experiments prob- 
ing finite-size scaling behavior in 4 He at the superfluid 
transition [17] . Moreover, accurate studies of the critical 
properties of trapped systems, to check universality and 
determine the critical exponents, require a robust control 
of the effects of the confining potential. In this respect, 



one may actually exploit trap-size scaling, using it to in- 
fer the critical exponents from the data, analogously to 
finite-size scaling techniques for the accurate determina- 
tion of the critical parameters, see, e.g., Ref. 

Helpful discussions with E. Arimondo, P. Calabrese 
and M. Oberthaler are gratefully acknowledged. 
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